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ABSTRACT 

Aims. We present the first results of an ambitious ground-based observation programme conducted on 1-4 meter class telescopes. Our 
sample consists of 1097 active and presumably young stars, all of them being optical (Tycho Catalogue) counterparts of ROSAT All- 
Sky Survey X-ray sources in the northern hemisphere. In this paper, we concentrate on the optically brightest (Vt < 9? 5) candidates 
(704 objects). We acquired high-resolution optical spectroscopy in the Ha and/or lithium spectral regions for 426 of such stars without 
relevant data in the literature. We describe the star sample and the observations and we start to discuss the physical properties of the 
investigated stars. 

Methods. We used a cross-correlation technique and other tools developed by us to derive accurate radial and rotational velocities 
and to perform an automatic spectral classification for both single stars and double-lined systems. The spectral subtraction technique 
was used to derive chromospheric activity levels and lithium abundances. We estimated the fraction of young single stars and multiple 
systems in stellar soft X-ray surveys and the contamination by more evolved systems, like RS CVn binaries. We classified stars on 
the basis of their lithium abundance and give a glimpse of their sky distribution. 

Results. The sample appears to be a mixture of quite young Pleiades-like and Hyades-like stars plus an older lithium-poor population 
probably born within the last 1-2 Gyr. Seven stars with a lithium abundance compatible with the age of IC 2602 (about 30 Myr) 
or even younger were detected as well, although two appear to be lithium-rich giants. The discovery of a large number of highly or 
moderately lithium-rich giants is another outcome of the present survey. 

Conclusions. The contamination of soft X-ray surveys by old systems in which the activity level is enhanced by tidal synchronisation 
is not negligible, especially for K-type stars. Five stars with lithium content close to the primordial abundance are probably associated 
with already known moving groups in the solar neighbourhood. Some of them are good post-T Tauri candidates according to their 
positions in the HR diagram. 

Key words, stars: activity - stars: fundamental parameters - stars: kinematics - stars: pre-main sequence - Binaries: spectroscopic - 
X-rays: stars 



1. Introduction 

Even if the global scenario of star formation history of the Milky 
Way i s fairly well known fro m its birth until the last few giga- 
years (lHavwood et al.lfl997h . surprisingly the recent local star 
formation rate is hardly constrained. Although newborn stars 
are usually confined to the locations of their parental clouds 
and share the same region of space (i.e. open clusters or asso- 
ciations), several processes tend to disperse them. Ejection of 
unused gas can reduce the mass of a cluster so much that it 
becomes gravitationally unbounded and breaks up. Moreover, 
occasional close encounters with giant molecular clouds along 
galactic rotation may even rapidly disrupt the cluster entirely. As 
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a result, most open clusters only survive a few hundred million 
years, although the cluster lifetime sensitively depends on the 
cluster's mass (for example M67 is a cluster with a solar age). 
Once mixed in the galactic plane environment, young stars are 
virtually indiscernible because neither their global photometric 
properties nor the presence of nearby gas can help to disentan- 
gle them from older stars. As a consequence, the recent (i.e. last 
giga-years) local star formation rate is poorly known because of 
the difficulties encountered in properly selecting young main se- 
quence (MS) late-type stars in the field from optical data alone. 
The same problem holds for the scale height of stars. As stars 
get older they tend to scatter higher above the galactic plane, a 
process known as disc heating. The thermalisation of the young 
stellar population is still an open question, and their scale height 
is generally assumed to be similar to the one of the molecular 
clouds from which they are born. 
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Fortunately, the situation has improved by making use of X- 
ray observations. The discovery that stars of almost all spectral 
types were X-ray emitters with luminosities in the range 10 26 - 
10 34 erg s 1 was on e of the major results of the Einstein mission 
dVaiana et al.l lT981). The current knowledge of the composition 
of the soft X-ray sky is mainly based either on compl ete optical 
identificatio ns of ROSAT All-Sky Survey ( RASS: IVoges et all 
19991 l2000l) sources in small-size regions dMotch et al.l Il997t 



Zickgraf et al] 1 19971 120051) or on the cross-correlation of the 



RASS with Tycho and Hipparcos catalo gues dESA|[l99 7). These 



last ones (namely RasTyc and RasHip; Guillou t et al. 19991) al- 



lowed study of the large-scale (all-sk y) distribution of X-ra y ac- 
tive stars in the solar neighbourhood (Guill out et al.ll 1 998ah and 
led to the disco very of the late-type stellar population of the 
Gould Belt (GB; lGuillout et al.lll998bl) . 

As chromospheric and coronal activity levels decrease 
with increas i ng age, they can serve as pro xy age indicators 
dHenrv et all 119961: ISterzik & Schmitj 1 19971) . During the last 
decade a fascinating picture of the recent star formation history 
in the solar neighbourhood has emerged. About fifty Myr ago, 
molecular clouds were forming stars near the present position 
of the Sun. Thirty Myr ago, the first generation of massive stars 
exploded as supernovae triggering star formation in an expand- 
ing ring-like structure and modelled the gas cavity in the solar 
neighbourhood (i.e. the local Bubble). Although the present rim 
coincides with most of the nearby OB associations, the distribu- 
tion of young stars as outlined by the RasTyc sample indicates 
that stellar formation also took place also inward the Belt, over a 
significant, yet poorly constrained, radial extent (the lack of ac- 
tive star formation within 50 pc of the Sun is believed to be due 
to the local Bubble). Although the massive-star content has been 
extensively studied, less is known about low-mass star formation 
in the Belt. Its age is uncertain by a factor of 2 because of the dis- 
crepancy found between the dynamical tim e scale (20-30 Myr; 
iMoreno et al]|l999t iPerrot & Grerrierll2003b and the s tellar age 
(30-60 Mvn lComeron & Torralll994tlTorra et al.ll2000b . 

It is known from statistical studies that open clusters account 
for only about 10 % of the total galacti c star f ormation. However, 
IZinneckerl d2008l) and iPiskunov et all (120081) have recently sug- 
gested that it may account for up to 50%, the remainder 
occurring in OB associations. Nevertheless, isolated very young 
stars have recently been detected. Although TW Hya displays 
all the characteristics of T Tauri stars (TTS), it has earned 
notoriety mostly on the basis of the absence of dark clouds in 
its vicinity. Different scenarios have been proposed as possible 
explanati ons for the presence of T TS far away from molecular 
clouds. Sterzik & Durise 3 d!995l) suggested that isolated pre 
main-sequence (PMS) stars have been ejected (run-away TTS) 
from their birthplace as a consequence of close encounters with 
other members of their parent cloud. Alternatively, they may 
have been form ed locally from small clouds, as suggested by 
Feigelson] d!996l) . As an isolated TTS, the origin of TW Hydra 
has long been a mystery f or as tronomers. However, s e arches 
by Ide laReza etaiT ( 1 1989b and iGregorio-Hetem et alj d 1992b 
have revealed four other TTS in the same region of the sky that 
shares the same space motion. At a distance of only » 50 pc 
and an age « 12 Myr, TW Hya is now recognised as one of the 
closest know n regions of recent star form ation, i.e. the TW Hya 
association (Zuckerman & Webb l2000h . Since its discovery, 
eight more young nearby associations have been identified so 
far dZuckerman et alJJ2006h . and eve n more are suspecte d on the 
basis of kinematics or other features dTorres et al.l 2008). 



Another issue regarding young stars is known as the post-T 
Tauri star (PTTS) problem. PTTSs cover the evolutionary period 
between 10 6 to 10 7 years (i.e. between the T Tauri phase and 
zero-age main sequence ZAMS age). They are expected to be 
abundant in the vicinity of star formation regions, although 
studies that concentrate on such regions often fail to find the 
number of PTTS that are expected fro m star formati on and 
evolution theories, as already note d by iHerbid dl978b thirty 
years ago. ISoderblom et al.l dl998bl) discussed possible expla- 
nations of the PTTS problem and in particular suggested that 
numerous PTTS may exist in isolated environments, but they 
are difficult to recognise because they mostly lack the obvious 
optical signatures of their T Tauri progenitors. 

Given their age and proximity, young, nearby associations 
and PTTS are important for understanding the evolution of 
circumstellar discs (in particular the transition phase of disc- 
dispersal) and the formation of planets, as well as for giving 
new insight into the process of planet migrations. RasTyc 
and RasHip samples may contain hundreds of such hidden 
young stars, and complementary optical data are absolutely 
necessary if the huge scientific potential of these samples is to 
be fully exploited. Moreover, detailed comparisons of radial 
velocity, proper motions, age, spectral type, and luminosity class 
distributions with pr edictions of stellar X-ray population model 
dGuillout et al.1 119961) are likely to establish new constraints on 
galactic evolution parameters such as the recent local stellar 
formation rate and the scale height of young stars. 

With this scientific background in mind, we started an opti- 
cal high-resolution spectroscopic campaign of a large sample of 
RasTyc / RasHip stars on 1- to 4-meter class telescopes. The 
sample is mostly composed of G and K late-type stars with 
Tycho B — V colour index in the range 0.6-1.33. The present pa- 
per is dedicated to the optically bright sample (composed of stars 
with Vt < 9?5) for which we already have a complete dataset. 
We describe the observations and discuss stellar properties as de- 
rived from high-resolution spectra. An outstanding result is the 
discovery of new members of moving groups and possible PTTS 
candidates, as well as the detection of several lithium-rich giants. 
We note that a comple mentary study in the southern hemisphere 
has been conducted bv lTorres et alj d2006l) . In the following we 
first introduce our stellar sample, the observations, and the re- 
duction methods (Sect.|2]i, then we present our determination of 
the stellar parameters (Sect. [3), the level of chromospheric activ- 
ity, and the age evaluation (Sect. [4]). In Section[5]we discuss the 
fraction of multiple systems, evolutionary status, age, and mov- 
ing group membership. We conclude and summarise the first re- 
sults from the ongoing study of these data (Sect. [6j. 



2. Sample selection and observing campaign 

As a starting point we used the RasTyc sample resulting from 
the cross-correlation of the RASS with the Tycho catalogue. To 
obtain statistically significant results and to fulfil our scientific 
goals, we selected about 1 100 RasTyc sources in the solar neigh- 
bourhood given the following constraints: 

- declination 5 > 0° ensuring that any target can be observed 
at low airmass from French and Italian national telescopes; 

- right ascension a e [15 H, 8 H] to cover both low and 
high galactic plane regions observable from the May-June to 
December-January periods decided as our observing slots in 
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Table 1. Observation program summary 



/■Si 




Fig. 1. Galactic distribution (Aitoff projection) of the stars se- 
lected for this spectroscopic survey (big circles) overlaid on 
the 8593 RasTyc stars with PSPC count-rate larger or equal to 
0.03 cts s (small dots). Only the optically bright sample (filled 
circles) is investigated in the present paper. The remainder (faint 
sample; open circles) will be studied in a subsequent paper. 



the framework of the Observatoire de Haute Provence (OHP 
- France) key programme; 

- 0.6 < (B - V)t < 1.3 range for which the lithium depletion 
behaviour is more effective. Indeed, it is well known that the 
strength of the Li 1/1670 7.8 line can be used to indicate youth 
(Martin &Clarej[i996) for stars cooler than about mid-G, 
although for M-type stars, the lithium is burned so rapidly 
in their deep convective envelopes that the Li 1/16707.8 line 
is only detectable for extremely young stars. The lithium 
plateau observed for F-type stars does not permit the use of 
this line as an age diagnostic; 

- ROSAT PSPC count rate greater t han or equal to 0.03 cts s" 1 
to avoid RASS scanning bias (see lGuillout et alJ [T999); 

- Vt < 10. n 5, since Tycho is largely incomplete above this 
limit. 

We show in Fig. Q] the galactic distribution of the make-up 
sample overlaid on the RasTyc all-sky picture, while Fig. [^dis- 
plays the selected stars magnitude and colour distributions. The 
median visual magnitude of the sample is ss 9^0. To gather as 
much as possible a-priori information on these potential targets 
and to minimise observing time, the data from the literature 
were checked against the SIMBAD database. Stars brighter than 
Vt ~ 6™0 mostly have over 50 bibliographic references with at 
least one spectroscopic investigation so were not scheduled for 
the observations. Those fainter than this limit were checked one 
by one and discarded only if relevant spectroscopic data already 
existed. Nevertheless, a few sta rs in common with the KPNO 
survey dStrassmeier et alj|2000l) have been observed to test the 
consistency of our results with previous works. 

The main observation programme was conducted mostly in 
service mode at the OHP 1 .52-m telescope from Nove mber 2001 
to Au gust 2005. We used the Aurelie spectrograph (Gille t etaT] 
Il994t) in both the Ha and lithium spectral regions with grating 7 
(18 000 grooves/mm), which allowed us to cover over 120 A at 
a resolution of R « 38 000. To derive reliable equivalent widths 
and good radial and rotational velocities, a signal-to-noise ratio 



Year 


Period 


Nights 


Telescope 


Spectra 


2000 


July - September 


22 


OHP 1.93-m 


60 


2001 


August 


12 


OHP 1.93-m 


50 


2001 


October - December 


29 


OHP 1.52-m 


290 


2002 


June - November 


53 


OHP 1.52-m 


396 


2003 


January 


23 


OHP 1.52-m 


207 


2004 


June - September 


29 


OHP 1.52-m 


231 


2005 


July - August 


20 


OHP 1.52-m 


126 


2007 


September - October 




OAC 0.91-m 


28 


2008 


March - July 




OAC0.91-m 


18 


Total 




188+ 




1 406 



Queue mode. 



Table 2. Instrumental setup summary 



Telescope 


Spectrograph 


Spectra 


Spectral region 


Resolution 


OHP 1.93-m 


Elodie 


echelle 


3900-6800 A 


42 000 


OHP 1.52-m 


Aurelie 


ID 


6495-6625 A 


38 000 


OHP 1.52-m 


Aurelie 


ID 


6650-6775 A 


38 000 


OAC 0.91-m 


FRESCO 


echelle 


4250-6850 A 


21000 



(S/N) of 70-100 per resolution element (depending on the star 
magnitude) was required. The Horizon 2000 CCD in bin mode 
permitted us to reach our requirements in the most efficient way. 
Another set of 1 10 stars observed in 2000 and 2001 (for a similar 
research project but over a smaller area) at the OHP 1.93-m tele- 
scope was merged with the previous sample. Spectra covering 
the 3900-6800 A range were acqu ired with the Elodie echelle 
spectrograph (Bar anne et al.l IT996b with a resolving power of 
R » 42 000. The equipment of both telescopes has about the 
same efficiency. The typical exposure times ranged from less 
than 10 minutes for a bright 7™0 star up to an hour for stars 
fainter than 9^0, for both instruments. Finally, 45 remaining stars 
brighter than Vt = 9^0, which could not be scheduled at the OHP 
because of time constraints or bad weather conditions, have been 
successfully observed at the M.G.Fracastoro station (Serra La 
Nave, Mt. Etna, 1750 m a.s.l.) of the Osservatorio Astrofisico di 
Catania (OAC - Italy). The OAC 0.91-m telescope was equipped 
with the FRESCO fiber-fed echelle spectrograph covering the 
spectral range 4250-6850 A with a resolution of R 21 000. 
Spectra of radial and rotational velocity standard stars, as well 
as bias, flat-field, and arc-lamp exposures, were acquired at least 
twice a night for calibration purposes. In total « 1400 spectra 
were acquired during the course of the ^190 nights devoted 
to the programme (Table Q]) with various instrumental setups 
(Table©. 

The ID Aurelie spectra were corrected for bias, flat-field, and 
cosmic rays hits by using standard MIDAS procedures while, 
the OAC data reduction was performed by using the ECHELLE 
task of IRAfQ. Elodie spectra were automatically reduced on- 
line during the observations with the standard pipeline reduc- 
tion package, which also includes the cross-correlation with a 
refer ence mask and the TGMET facility for spectral classifica- 
tion dSoubiran et al.l 1998b . Tell uric lines were sub tracted follow- 
ing the procedure described bv lFrasca et alj d2000l) . and spectra 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association for research 
in Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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Fig. 2. Tycho Vt magnitude (left panel) and (B - V)j colour (right panel) distributions of the stars in our sample. Stars with 
already available relevant measurements from the literature (green right-stripped area) have been systematically discarded from the 
observation programme. Spectroscopic observations of the optically bright sample have been conducted mostly on the OHP 1.52-m 
telescope (grey area), while a significant part of the optically faint sample has been scheduled on the larger OHP 1.93-m and TNG 
3.58-m telescopes (red left-stripped area). A small number of stars (white area ; Vt < 9?5) have also been scheduled on the OAC 
0.9 1-m telescope. 



were then normalised to the continuum through a polynomial 
fit avoiding spectral regions heavily affected by blends of ab- 
sorption lines, molecular bands, and the Ha wings. Spectra were 
finally corrected for the Earth velocity and reduced to the helio- 
centric rest frame. 

3. Determination of astrophysical parameters 

3.1. Multiplicity, radial and projected rotational velocity 

For each star observed with Aurelie, we computed the cross- 
correlation function (CCF) independently for both the Ha and 
lithium spectral regions by cross-correlating the observed spec- 
tra with a synthetic mask. Since the CCF improves when spectra 
with a very similar shape are correlated, we used a grid of 9 dif- 
ferent synthetic spectra (corresponding to F0, F5, GO, G5, K0, 
K5, MO, and M5 spectral types) as CCF masks, and selected the 
one that is more similar to the spectral type (or the B - V colour 
index) of the target. The CCF is a powerful tool for revealing 
single, double, and multiple systems. Significant minima in the 
CCF were searched for automatically and single, double, or mul- 
tiple functions (i.e. up to 4 Gaussian or rotational profiles) were 
fitted to the CCF dips (see Fig . [3j . The standard deviation of the 
CCF values far from the central dip(s) (<x) was evaluated for each 
spectrum and all the minima deeper than 3<x were kept as signif- 
icant. This allowed us to sort out single stars (S), double-lined 
binaries (SB2), triple systems (SB3), and multiple stars (M). 

In principle, the minimum of the CCF readily gives the radial 
velocity (RV) of the star. However, the CCF profile may appear 
asymmetric or distorted because of cool star-spots on the stellar 
photosphere. Therefore, we fitted the entire CCF profile to obtain 
a more accurate RV determination. 



For most of our spectra, the RV uncertainty is dominated 
by the accuracy of the wavelength calibration. The effects of 
spectrograph drift during the night were minimised by acquiring 
arc lamp exposures at least twice per night. An overall RV 
precision of + 1.5 kms -1 was estimated by comparing our RV 
measurements with those from the literature for a set of standard 
stars (see Fig. |4]i although the precision can reach 0.5 km s~ l for 
high S/N spectra. 

We also took advantage of each star normally being observed 
a few nights apart, in the Ha and Li i regions, to gather spectro- 
scopically unresolved binaries. For stars appearing as singles, 
i.e. with single-lined spectra without any detectable RV varia- 
tion on a timescale of a few days, we averaged the two RV val- 
ues. On the other hand, significant RV differences between the 
He and lithium spectra were interpreted as real variations. Thus, 
the stars with a single dip in the CCF and a variable RV were 
classified as single-lined binaries (SB1; see Fig- [5j whenever the 
differences in RV measurements were greater than the sum of 
their uncertainties. 

The CCF also permits the projected rotational velocity 
(v sin i) to be determined. The shape of the CCF is sensitive to 
the rotation speed. At low v sin i, i.e. when the spectral lines are 
mainly affected by the Doppler and microturbulence broaden- 
ings, a Gaussian function generally gives a very good approx- 
imation of the CCF profile and we used a calibration relation 
between the width of the fitted Gaussian, <x(CCF), and v sin i to 
determ ine the rotational ve locity, following a method similar to 
that of lOueloz et alj (119981) . At high v sin i values, the shape of 
the CCF is far from Gaussian and becomes more and more like a 
rotational profile. In such cases, a Gaussian function would over- 
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Fig. 3. Cross-correlation function (asterisks) and Gaussian fitting (red solid line) for the single star HD 6026 (RasTyc 0102+6236; 
left panel), the SB2 binary BD+67 34 (RasTyc 0025+6848; middle panel), and the triple system BD+82 622 (RasTyc 2034+8253; 
right panel). Lower panels show residuals. 




RV (km/s) - Literature 



Fig. 4. Comparison of our radial velocities measurements with 
those quoted in the literature for the standard stars observed 
during our runs. Residuals from the one-to-one relation are 
shown on the bottom panel. The short-dashed lines delimit a 
±1.5 km s region around the one-to-one relation (continuous 
line in the upper box). 



estimate the rotational velocity, so we preferred to compare the 
observed spectra with a grid of synthetic templates broadened 
by the convolution with a rotation profile and to retain the v sin i 
giving the lowest residuals in the difference. 



< / 



Fig. 5. Radial velocities measured from the lithium (RV/,,) ver- 
sus Ha (KVh„) spectral regions for stars classified as single from 
the CCF profile analysis. The dashed line shows the one-to- 
one relation. Stars with inconsistent RVl,-RV#„ measurements 
(filled red circles) are interpreted as single-lined spectral bina- 
ries (SB 1). 



The error on the projected rotational velocity was computed 
for each star and depends critically on S/N ratio, as well as 
on the shape of the cr(CCF) versus vsim relation. At low 
vsinz (< 25 km s _1 ), the relation can be approximated with a 
square-root function and a small error on the <x(CCF) propagate 
on a significant vsim relative error (50 %). On the other hand, 
at high v sin i, the relation is rather linear and the relative error 
is lower (10 %). In any case, our spectral resolution does not 
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permit to derive v sin i values lower than about 5 km s _I . Finally, 
one must keep in mind that unresolved lines from a secondary 
companion can mimic a higher projected rotational velocity 
(and can bias the estimation of the astrophysical parameters). 

For the 45 stars observed with FRESCO at OAC, we adopted 
spectra of F-G-K radial velocity standard stars acquired during 
the observing runs as RV templates for the cross-correlation. 
We used all the 20 echelle orders, excluding the Ha and Nai 
D2 lines from the analysis, because of their strong wings affect- 
ing the CCF shape and due to the contamination of these lines 
by chromospheric emission. The radial velocity measurements, 
performed with the IRAF task Fxcor, were obtained by averag- 
ing RV data from all echelle orders with the usual instrumental 
weight Wi = err 2 . The cr ; values were computed by Fxcor ac- 
cording to th e fitted peak height an d the antisymmetric noise as 
described by lTonry & Davisl (Il979l) . The standard errors of the 
weighted means were computed on the basis of the errors en 
in the RV values for eac h order according to the usual formula 
(see, e.g.. lToppinglll957h . For the v sin i determination, we used 
the same method as for the Aurelie spectra but the spectra of 
slowly-rotating standard stars were used as templates. 

3.2. Effective temperature, gravity, and metallicity 

In the case of single stars or SB1 systems (i.e. only one signif- 
icant minimu m in the CCF; see Fig. [3] left panel), we run the 
ROTFIT code jFrasca et al.| [2003). which finds the spectral type 
of the target searching into a library of standard-star spectra, 
the spectrum that gives the best match (minimum of the resid- 
uals) with the target one (see, for instance, online Fig. ITTb. after 
rotational broadening (hence giving another independent mea- 
sure of the star's vsirn"). As a standard-star library, we used a 
set of 1 85 spectra of stars w ith known astrophysical parameters 
( Prugniel & Soubiran 200 1 ) retrieved from the ELODIE Archive 
dMoultaka et alj 120041) . which are well-distributed in effective 
temperature and gravity and in a suitable range of metallicities. 

For double stars, we used anoth er method (implem ented as 
an IDL code), namely COMP02 dFrasca et alj 120061) , which 
searches for the best combination of two standard-star spectra 
able to reproduce the observed spectrum of the SB2 system. 
We give, as input parameters, the radial velocity and v sin i of 
the two components, which were derived from the CCF anal- 
ysis. COMP02 then finds, for the selected spectral region, the 
spectral types and fractional flux contributions that better repro- 
duce the observed spectrum, i.e. which minimise the residuals in 
the collection of difference (observed - synthetic) spectra. For 
this task we used a smaller set of 87 reference spectra retrieved 
from the ELODIE Archive that are representative of stars with a 
nearly-solar metallicity, spectral types from early-F to early-M, 
and luminosity classes V, IV, and III. An example of the results 
provided by COMP02 is shown in Fig. [6] 

For multiple systems, we made no attempt to derive as- 
trophysical parameters, and these systems will be analysed in 
detail in a future dedicated work where a new COMPO code, 
designed for multiple systems, will be used. 

The accuracy of the stellar parameters (T e ff, logg, [Fe/H]) 
depends not only on the S/N of the spectra but also on the homo- 
geneous coverage of the grid of parameters by reference stars. 
To check the reliability of our determinations, we compared our 
results with those reported in the literature for six binaries clas- 
sified with COMP02. The determination of stellar parameters 
for the components of SB2 systems is a challenging task for any 
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Fig. 6. Observed Aurelie spectrum (black thick line) of the SB2 
system BD+15 2915 (RasTyc 1547+1509, K2V+K3V) in the 
Ha (upper panel) and lithium (lower panel) spectral regions, to- 
gether with the synthetic spectrum built up with the weighted 
sum of two Elodie reference spectra broadened at the v sin i of 
the system's components and Doppler-shifted according to its 
RV (thin red line). The Ha cores of the primary (P) and sec- 
ondary (S) components are marked in the upper panel, while the 
Cai /16717 lines of the two components are labelled in the same 
way in the lower panel. 



automatic spectral classification code, due to the blending of the 
spectral lines of the two stars. As COMP02 is the equivalent of 
the ROTFIT code adapted to binary systems, we consider the ac- 
curacy on the parameters of SB2 components as an upper limit 
on the uncertainties for single stars. As apparent from Table [3] 
our spectral type classification agrees with that from the litera- 
ture within 1-2 subclasses for nearly all systems. 

As a second test, we selected stars classified by ROTFIT 
for which the trigonometric distance has an uncertainty below 
30% (D30) and plotted them on an HR diagram, as illustrated 
in Fig. [7] left panel. It readily shows that none of the stars 
lying on the MS was misclassified as a giant and that, among 
the 50 stars with B - V > 0.8, only 6 (i.e. 12 %) giants 
or subgiants (according to their HR diagram location) were 
misclassified as dwarfs. On the other hand, stars leaving the 
upper main sequence, but which have not yet reached the giant 
clump, are systematically classified wrongly as dwarfs by our 
automatic procedure. Quantitatively, 22 among the 111 stars 
with B - V < 0.8, i.e. = 20%, are misclassified. This is likely 
the result of the lack of similar stars in the reference library (due 
to the very short lifetime of this evolutionary phase), but not 
because of the formal algorithm of our code. 
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Table 3. Stellar parameters for 7 binaries with literature values. 



RasTyc Name 


Name 


Sp.Type" 1 " 


Sp.Type L " 


T Our 
J rff 




RasTyc 0702-0515 


AR Mon 


Kl III/G2 V 


.../... 




4395 / 5350 


4500 / 5200 1 


RasTyc 0921-0640 


NYHya 


G2 V / G2 V 


.../... 




5450 / 5450 


5490 /... 2 


RasTyc 1310+3556 


RSCvn 


F7 V / G8 V 


F6 IV /G8 


IV 3 


6120/5420 


.../... 


RasTyc 1601+5120 


EVDra 


G7 V/Kl V 


G8 V/Kl 


V 4 


5400 / 5370 


.../... 


RasTyc 1825+1817 


AW Her 


G2 V / K0 IV 


G2 IV / K2 


IV 5 


5780/4870 


.../... 


RasTyc 1912+4619 


FLLyr 


G2 V/Kl V 


F8 / G8 ( 




5730 / 5080 


6080 / 5390 7 



Re ferences. (lllWilliamon et alj d2005|): (2)IClausen et alJd200lh : Gl lStrassmeier & Fekell d!990h : (4) lHenrv et all dl995h : (5) lrmbertl 12002): 
(6) iPopper et alJdl986l) : (7)lMolenda-Zakow icz et alJd2007h l 



Table 4. Consistency of the stellar parameters derived from the Aurelie (left) and Elodie (right) spectrographs. 



RasTyc Name 


RV (km s-') 


v sin i 


(km s _1 ) 




(K) 


l°g£ 


[Fe/H] 


RasTyc 2100+4530 


-23.3 -23.1 


6 


6 


5720 


5730 


4.23 4.42 


-0.01 0.01 


RasTyc 2114+3941 


-10.8 -11.2 


24 


27 


5730 


5730 


4.14 4.28 


-0.23 -0.53 


RasTyc 2118+2613 


-16.3 -16.8 


7 


6 


5230 


5330 


4.31 4.52 


-0.09 -0.10 


RasTyc 2132+3605 


-6.2 -6.8 


5 


4 


5140 


5000 


3.20 2.84 


-0.03 -0.05 


RasTyc 2137+1946 


-5.2 -6.0 


4 


3 


4900 


4870 


3.03 2.85 


0.04 -0.01 



Thus, we consider our spectral classification affected by 
an uncertainty that has a typical value of one spectral sub- 
class (AT e ff « 150 K) but can reach two spectral subclasses 
(Ar e ff * 300 K) in the less favourable cases (especially for 
SB2 systems). We also stress that = 20% of our stars bluer 
than B — V — 0.8 can suffer luminosity misclassification; i.e., 
they could be stars on their way towards the giant branch that 
have a lower gravity compared to the MS stars with the same 
temperature. 

The data acquired from the various telescopes/spectrographs 
were treated in slightly different ways. That is why a few sin- 
gle stars were observed with different instruments for check- 
ing the consistency of our results. As apparent from Table |H 
which compares parameters derived from the T1.52+Aurelie 
and Tl. 93+Elodie spectrographs, there is no discrepancy or 
systematic effects and the agreement is very good on aver- 
age. Consistency of measurements fro m OHP and OAC_ spec- 
tra has already been demonstrated in Frascaetal. (2006) and 
iKlutsch et al.lO008h . 



4. Chromospheric emission and lithium abundance 

4. 1 . Chromospheric emission 

In addition to the fundamental parameters (effective tempera- 
ture, gravity, metallicity, and v sin /), the ROTFIT and COMP02 
codes allowed us to derive stellar activity and age. We es- 
timated the chromospheric emiss ion level with the "spectral 
subtraction" techn i que ( s ee, e.g., lHerbidfT98l lBardenlfl98l 
IFrasca & Catalanol Il994t lMontesetal.ll 1995b . which is based 
on the subtraction of a reference "non-active" spectrum built 
up with synthetic model-atmosphere spectra or with observed 
spectra of slowly-rotating stars of the same spectral type with 
a negligible level of chromospheric activity. The net equivalent 
width of the Ha line, W^™, was measured in the spectrum ob- 
tained after subtracting the non-active template by integrating 
the residual Ha emission profile. Although it is not easy to get 
truly chromospherically non-active stars, especially for the K- 
M spectral types, using the spectra exhibiting the minimum Ha 
residual flux as "non-active" templates allowed us to pick up 



stars with a moderate activity and to evaluate its level. The Ha 
luminosity, which is a parameter more indicative of the chromo- 
spheric activity level, was calculated from the net Ha equivalent 
width, the distance d, and the continuum Earth flux at the Ha 
wavelength, /6563, according to the following equation: 

L Ho = 4Kd 2 f 65 ^WZ 

= 4 Kd 2 ^10 ( - QAV °- &A5 "WZ, (1) 
-F5556 

where io(-°- 4y °- 8 - 451 > is the Earth flux at 5556 A of a star with 
de-reddened mag nitude V (iGrayl 19921) and the continuum flux- 
ratio was evaluated f rom NextGen synthetic low -resolution 
spectra ' lHauschildtetalJ Q999) at the temperature found with 
the ROTFIT code. 

4.2. Lithium abundance 

Lithium is strongly depleted when mixing mechanisms pull it 
deeply into the convective layers of late-type stars. A strong 
Li 1/16707.8 photospheric absorpti on line is thus generally con- 
sidered as a youth indicator (e.g.. ISoderblom et alJ [T998a). We 
measured the equivalent width of the Li 1/16707.8 line on the 
residual spectrum, once it hs been checked that this feature is 
not present in the spectrum of the reference star. The subtrac- 
tion technique allowed us to remove the contribution of the 
small nearby Fe 1/16707.44 line, whose intensity is expected to 
be nearly the same in the reference spectrum that best represents 
the target. We estimate the lower limit for positive detection of 
the Li 1/16707.8 line to be approximately 10 mA, although it 
strongly depends on the S/N of the spectrum. The lithium abun- 
dance was finally deduced from the Li 1/16707.8 line equivalent 
width (EW(Li)) and effect ive temperature by linear int erpolation 
of the values tabulated by Pavlenko & Magazzu (119961) . 

5. Discussion 

Although a detailed multidimensional analysis of all parameters 
derived from our spectra is beyond the scope of this paper, in the 
following subsections we briefly discuss the fraction of multiple 
systems (Sect. 15. It , as well as the evolutionary status (Sect. 15 - 2b 
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Fig. 7. Le/f panel: stars with trigonometric distance 7r known better than 30% (D30) and classified as class V (white square symbols) 
or class III, IV (black triangle symbols) superposed on an observational HR diagram of the * 3 400 RasHip stars (dot symbols ; 
grey scale code their density gradient) with accurate distance (cr n /n < 0.1) and colour (ctb-v < 0.025). Right panel : Zoom 
on the late-type region of the HR diagram showing the main sequence and giant clump regions together with D30 stars classified 
as PMS-like (filled circle ; #i show cross-match with Table |5), Pleiades-like (filled square) and UMa-like (open square) based on 
their lithium abundance. Siess et al. (1997) PMS (from 0.8 to 1.2 M ) and post MS (from 1.2 to 3.0M Q ) evolutionary tracks and 
isochrones (1,3, 10, 30, and 50 Myr from top to bottom) have been overlaid. Zero age main sequence (bold solid line) and terminal 
age main sequence (dashed bold line) are shown as well. The 2 PMS-like stars are discussed in detail in Sect. [~ 



and age classification (Sect. \5.3\ of our star sample. Finally we 
investigate in detail the lithium-rich stars (Sect. |5.4] and Sect. 15. 5t 
of our sample, among which five appear as possible members of 
young stellar kinematical groups in the solar neighbourhood. 

5.1. Binaries and multiple systems 

Bran dner et al.l dl996h have pointed out that ROSAT unresolved 
binaries are X-ray sources statistically brighter than single 
stars. Thus, X-ray flux-limited surveys suffer from a detection 
bias towards binarity. Indeed, a large number of spectroscopic 
binaries (SB) have been discovered in this sample. More quan- 
titatively, if restricting ourselves to stars observed with Aurelie 
in both spectral regions, our sample contains 210 single stars, 
47 single-lined spectroscopic binaries (SB1), 38 double-lined 
spectroscopic binaries (SB2), and 9 bona-fide triple systems, 
among which 2 are potential multiple (more than 3 components) 
systems. Most of them are new detections (according to the 
SIMBAD database). Thus, in our sample we have percentages 
of 15.5%, 12.5%, and 3.0% for SB1, SB2, and multiple systems, 



respectively. Altogether, non-single stars account for 3 1 % of the 
sample. 

Six binaries ( 3 SB1 and 3 SB2;lFrasca et al.ll2006l) as well as 
3 triples systems (iKlutsch et al.ll2.Q08h have already been studied 
in great detail. The high-quality RV curves have been solved for 
these systems, and the orbital and physical parameters obtained. 
Although no significant conclusions can be drawn with respect 
to evolutionary status of SB, we are continuing to monitor a rep- 
resentative sample of RasTyc SB. 

5.2. Evolutionary status 

Stellar soft X-ray surveys have long been assumed to be mostly 
composed of young active late-type stars with a possible con- 
tamination by older active binaries (i.e. RS CVn systems; 
ISciortino etalH l995). However, because of the lack of knowl- 
edge on this population (density and scale height), the compari- 
son of stellar X-ray population models with optically-identified 
X-ray sources generally assumed that the contamination by 
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Fig. 8. Li 1/16707.8 line equivalent widths versus B-V colour 
index for stars of our sample. The lines represent the boundaries 
for IC 2602 (orange dash), Pleiades (red solid), and Hyades (blue 
dash-dot) clusters. The black long-dashed line shows our lower 
limit for positive detection of the Li 1/16707.8 line. Seven stars 
(red dots next to #i symbols; see Table|6]for cross-match) show 
lithium abundance in excess of the lithium-richest Pleiades stars 
and 2 are located above IC 2602 lithium upper envelope. 



evolved systems is at most only marginal. Depending on its 
strength, the impact of the contamination by older systems that 
mimic young stars can lead to significant error o n the derived re- 
cent local star formation rate (Afferet al. 2008). Assuming that 
the RasTyc sample is representative of X-ray unbiased and mag- 
nitude limited surveys, we derived an overall contamination by 
evolved stars (likely RS CVn systems) of about 35 % with a peak 
of 60 % in the K-type stars region (0.8 < B - V < 1.1), whereas 
this fraction decreases to only 20 % and 10 % for hotter (F-G) 
and cooler (M) stars, respectively. 



5.3. Age 

The lithium abundance to indicate the empirical age is now 
commonly accepted for stars cooler than B - V - 07 1 6. It 
is possible to define with high confidence in a colour versus 
Li 1/16707.8 line equivalent width diagram the loci occupied by 
PMS-like, Pleiades-like, UMa-like, Hyades-like, and Old (i.e. 
age greater than one Gyr) stars. Such a diagram is presented 
in Fig. [8] for only our single an d SB1 stars , together with 
the IC 2602 dMontes e tail 12001 al) . Pl eiades (ISoderblom et all 
119931: lNeuhauserf[l997l) . and Hyades (ISoderblom etalj 1 19901) 
envelopes. Some of our stars exhibit a lithium content larger 
than or compatible with that of IC 2602, indicating PMS or 
PTTS candidates (see discussion in Sect. 15.51 ). Other stars 
occupy a domain similar to that of the Pleiades, UMa, or 
Hyades, so they likely have similar ages. There are also stars 
with no photospheric Li 1/16707.8 line detected or with lithium 
content below the Hyades lower envelope, which are likely to 
be slightly older than one giga-year according to current stellar 



X-ray population models. 

The fractions of PMS-like, Pleiades-like, UMa-like, Hyades- 
like, and Old stars are 2.6%, 7.0%, 39.1%, 19.2%, and 32.1%, 
respectively. Interestingly, the fraction of PMS-like stars 
derived fro m the only lithium criterion is exactly the one 
predicted by iGuillout et al.l (|T999) from a statistical analysis of 
RasTyc/RasHip samples. Altogether, the very young stars (i.e. 
PMS-like plus Pleiades-like and UMa-like) contribute about 
half of the sample. If we consider all the young stars (i.e. very 
young plus Hyades-like, namely stars younger than about one 
giga-year), this fraction increases to about 70%. Outside star 
formation sites (clusters or associations), young stars generally 
contribute only marginally (< 10% assuming a constant star 
formation rate) to optical star counts, in contrast to the high 
fraction found by us. Although lithium is depleted very rapidly 
in the convective envelopes of cool stars, it is not surprising to 
detect a significant fraction of lithium-rich stars in our sample, 
which is heavily biased towards young stars thanks to X-ray 
selection. 

A crude analysis of the galactic latitude distribution of 
young and old populations of RasTyc stars does not reveal any 
enhancement of young star number towards the galactic plane. 
This is actually not surprising because the majority of our stars 
(70%) are located within 100 pc of the Sun, a radius that is 
smaller than or comparable to the scale height of young stars 
estimated in between 70 pc (assu ming the scale height of molec- 
ular clouds : ISanders et al.lll984|) and 140 pc (as computed in the 
Besan con synthesis model ; Robin & Crezd [T986T: iRobin etail 
2003). The comparison of our sample with the pred ictions of a 
stellar X-ray population model (Guillo ut et al.l [l996) will be the 
scope of a dedicated paper. 



5.4. Lithium-rich giants 

In the right panel of Fig. [7] we plot a zoom on the late-type 
region of the HR diagram showing the main sequence and giant 
clump regions, together with D30 stars classified as PMS-like, 
Pleiades-like, and UMa-like based on their lithium abundance. 
It shows that, unforeseeasly, a significant fraction (23 out of 
50, i.e. 46 %) of the evolved stars are also lithium-rich. Most 
of them are classified as UMa-like but a few of them display an 
even larger lithium content and are classified as Pleiades-like 
or PMS-like, accordingly. These stars have early-F or A-type 
progenitors on the main sequence (M 1.2-3.5 Mq), which 
are known to be very fast rotators but which generally do not 
display any Li 1/16707.8 absorptio n line feature, although some 
Ap stars have lithium-rich spots dDrake et al.ll2004l 120051) . 



(Ide La Reza et al.| 1 19961 ll997allbL l2000t iDrake et al.l 120021: 
proved that all K-giant stars, with masses 



In a s eries of works, de la Reza and co-workers 

etal.) 

Ide La Rezal 12006b 
between 1 and 2.5 solar masses, become lithium-rich during the 
red giant branch stage. For these authors, these stars are normal 
giants going through a short period of lithium enrichment. 
This episode, which could be cyclic, would take origin in the 
changes of the stellar structure of the star. As these stars evolve 
towards the blue part of the clump (as helium-burning giants), 
they develop outer convection zones and consequently a high 
magnetic activity level, because of the convection and fast 
rotation (Fekel & Balachandran 1993), which yields the X-ray 
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Table 5. PMS-like candidates based on the lithium criteria only. 



# Name 


RasTyc Name 


a (2000) 


6 (2000) 


PSPC 


Sp. T. 




l°g£ 


[Fe/H] 


log(n L ,,) 


log(Hor) 






h m s 


o ' " 


(ct/s) 


L. cl. 


(K) 










#1 BD+45 598 


RasTyc 0221+4600 


02 21 13.0 


+46 00 7.6 


0.35 


Kl V 


5149 


4.04 


-0.09 


3.20 


28.97 


#2 BD+29 525 


RasTyc 0307+3020 


03 07 59.2 


+30 20 26.7 


0.45 


G2 V 


5827 


4.31 


-0.03 


3.24 




#3 HD 275148* 


RasTyc 0319+4212 


03 19 52.8 


+42 12 31.1 


0.03 


G7III 


5158 


3.39 


-0.00 


3.66 


30.46 


#4 HD 22179 


RasTyc 0335+3113 


03 35 29.9 


+31 13 37.8 


0.26 


G3 V 


5705 


4.32 


-0.06 


3.13 


29.37 


#5 HD 23524 


RasTyc 0348+5202 


03 48 23.0 


+52 02 16.9 


0.86 


G9 V 


5213 


4.22 


-0.06 


2.94 


29.33 


#6 HD 170527* 


RasTyc 1825+6450 


18 25 10.1 


+64 50 17.8 


0.23 


G4III 


5118 


3.60 


-0.09 


3.69 


31.10 


#7 BD+44 3670 


RasTyc 2100+4530 


21 00 47.0 


+45 30 10.4 


0.35 


G2 V 


5719 


4.23 


-0.01 


3.21 


29.26 



* were later discarded as PMS candidates because of spectral classification and kinematic criteria inconsistancies. 



emission. 

The mechanism leading to simultaneously high rotation, 
high lithium abundance, and enhanced act ivity is not yet un- 
dersto od. A first scenario was proposed by ICameron & Fowlerl 
dl97ll) for the asymptotic giant branch evolution phase: the 
convection zone can go so deep inside the stellar interior 
that it carries material cont aining 7 Be to upper c ooler layers 
where it decays into 7 Li. iFekel & Balachandra ^ (119931) have 
proposed a scenario in which internal convective motions 
penetrate into the rapidly rotating core dredging up to the 
surface high angular mo mentum material and lit hium recently 
synthesised. In contrast. Ide La Reza et alJ (|l996) suggests that 
the enrichment of lithium is an obliged and rapid phase for 
all the K-type g iants, produced by a mech anism similar to the 
one invoked by ICameron & Fowlerl d 197 lb . which leads to the 
formation of a shell of gas and dust. Other studies suggest 
that the lithium enrichment could originate in th e ingestion of 
subst e llar bodies like planets or brown dwarfs dSiess & Livid 
119991 iDenissenkov & Weissl [2000). At this time, the Lii line 
should appear in their spectra. However, the mechanism of 
lithium formation loses its efficiency after some 10 5 years so 
that thereafter the lithium is progressively destroyed in the 
deeper layers where it is brought by the convecting mixing. 
The observation of a strong Li 1/16707.8 absorption line would 
indicate that these stars have just reached this evolutionary 
phase without having enough time to deplete the photospheric 
lithium. 

Last but not the least could be lithium preservation in 
atmospheric stellar layers after some process like those acting in 
Ap stars (high and closed magnetic fields, diffusion, etc.) when 
the stars evolve out of the main sequence; but this scenario 
should suggest that lithium-rich K-giants may have spotted 
lithium-rich Ap progenitors, a new scenario that needs more 
detailed analysis in the future. 

The discovery of this large number of lithium-rich gi- 
ants adds more than 40 % to the kno wn ones (40 quote d by 
ldeLaRezaeTaT]|1997bl 15 quoted by iPerevra etafl 120061 and 
some oth er dispersed in the lit e rature ). Our results reinforce 
those of IFekel & Balachandraril d 1994b . who found that 58 % 
of single active giants have high lithium abundance (36 % for 
the giants in binary systems) and are also in line with a 50 % 
of lithium-r i ch sta rs among rapidly-rotating K-giants found by 
iDrake etail ([2002). The emerging picture is that a high or mod- 
erate atmospheric lithium abundance, coronal activity, and a 
rather rapid rotation are characteristic of an evolutionary phase 
that should be relatively short, given the small fraction (1-2%) 



of lithium- rich stars among the whole sample of G- and K-type 
giants (see lCharbonnel & Balachandran 2003). 

5.5. Moving groups candidates 

The highest EW(Li) values measured in our sample reach 
over 350 mA for 2 stars, namely RasTyc 1825+6450 
(#6) and RasTyc 0319+4212 (#3). Together with 5 other 
stars (RasTyc 0221+4600, #1; RasTyc 0307+3020, #2; 
RasTyc 0335+3113, #4; RasTyc 0348+5202, #5 and 
RasTyc 2100+4530, #7) showing EW(Li) compatible with 
that of IC 2602, they have a lithium content close to the 
primordial abundance log(n L(/ ) -3.2 (see Table [5} and can be 
considered as trustworthy PMS candidates. The observed spec- 
tra of our 7 candidates are presented in Fig. QT|and Fig- El It 
readily shows that, although all sources (with the only exception 
of RasTyc 0319+4212) display significantly filled-in Ha line 
profiles typical of stellar activity, none of them exhibits strong 
emission features above the continuum that are distinctive of 
classical T Tauri stars. On the other hand, no spectroscopic 
criterion has yet been clearly ident ified to disenta ngle weak-line 
T Tauri stars (WTTS) from PTTS. lMartfnld 1998b has tentatively 
established a divided line in the EW(Li) - diagram for 
distinguishing between WTTS and PTTS but unfortunately our 
candidates fall in the region of confusion where uncertainties 
prevent us from distinguishing them. 

Although RasTyc 1825+6450 and RasTyc 0319+4212 are 
the lithium-richest stars of our sample, they were classified 
as giants by our automatic spectral classification code (see 
Sect. |3.21 l. A strong lithium absorption in the atmosphere of 
giants stars is not an obvious indicator of youth, as we shortly 
discussed in Sect. 15.31 However, very young PMS stars that are 
still in the radiative part of their pre-main sequence evolutionary 
track can mimic evolved (subgiant or giant) stars with typical 
gravities in the range log g = 3.0-3.5, so our code can find 
a better match with a subgiant or giant template. In addition, 
their position in the HR diagram lies in the region occupied by 
subgiants of higher mass. Thus, we cannot reject them from the 
list of PMS star candidates. 

The positions (galactic coordinates) of candidates #1, #2, 
#3, #4, and #5 are overplotted in F ig. |9]on a map of dark clouds 
published by Dob ashi et alJ d2005b . The closest prominent star 
forming region is the Taurus-Auriga complex located more than 
10 degrees east of star #4 outside the portion of the displayed 
map. Although this candidate is projected directly in front of 
the Perseus molecular c loud, in the vicinity of the Per OB2 
association dBlaauwlll952b and IC 348 young star cluster, its 
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Fig. 9. Position (galactic coordinates) of PMS-like candidates overplotted on a map of dark clouds enhancing interstellar material 
(see #i next to cross symbols for cross match with Table |6]l 



distance is hardly compatible with these background objects, 
even assuming that the star is overluminous with respect to 
MS calibration because of its PMS nature. Thus, we consider 
this star as unrelated to the Perseus star-forming region. Star 
#5 is projected in front of unrelated background dark clouds 
on the galactic plane, while stars #1, #2, and #3 are located 
below the galactic plane in regions free from dense interstellar 
matter. Star #7 (not shown) is located just north of the unrelated 
North America nebula, while star #6 (not shown) lies above the 
galactic plane in the Draco constellation. Since none of them 
is located towards any close (i.e. within 200 pc) prominent 
star-forming region, they appear to be likely members of young 
stellar kinematic groups (SKG) in the solar neighbourhood. 

To assess possible membership as already known SKG, we 
computed the space velocity components (U, V, W; see Table [6) 
for these seven stars in a lefthanded coordinate system using 
Hipparcos/Tycho parallaxes (#5 and #6), otherwise photomet- 
ric d i stanc e derived from absolute magnit ude calibration from 
ICoxl d2000l) and Tycho2 (iFtog etal.1 12000) proper motions to- 
gether with our measured radial velocities. The U-V and U-W 
kinematic diagrams are shown in Fig. \W\ for stars listed in 
Table [5] We also show the location of stars satisfying one of the 
two Egge n criteria for t he ma jor stellar kinematics groups dis- 
cussed in Mo ntes et alj d2001bl) . namely the IC 2391 superclus- 
ter (~ 50 Myr), the Pleiades SKG (~ 100 Myr), the Castor SKG 



(~ 200 Myr), the Ursa Major (UMa) group (~ 300 Myr), and the 
Hyades supercluster (~ 600 Myr). The locus of the young-disc 
(YD; age < 2 Gyr) and the Sun position are shown as well. 

Figure[10]clearly shows that five stars have U, V, and W com- 
patible with the IC 2391 and/or Pleiades and/or Castor kinemat- 
ics, which are the youngest moving groups considered in this 
study. RasTyc 0319+4212 (#3) is at most marginally compati- 
ble with the Pleiades SKG kinematics while RasTyc 1825+6450 
(#6) is located far outside of the YD locus in the U-V and 
V-W planes and displays kinematics inconsistent with any 
SKG. Amazingly, these two stars are the ones classified as gi- 
ants by our automatic spectral classification programme. 

The membership in SK G can be analy sed more objectively 
using probabilistic methods (Klutsch 2008; Klutsch et al. 2009, 
in preparation). We computed the probability that a star with 
heliocentric velocities U, V, and W has a galactic motion com- 
patible with one of the five aforementioned comoving groups 
of stars. Results are presented in Table [6] We accepted a star 
as a possible new member of an SKG and flagged it with a 
Y in Table [6] if the probability (bracketed near the acceptance 
flag) exceeded 10%. According to this criterion, we find 5 stars 
displaying galactic motion compatible with the IC 2391 and/or 
Pleiades and/or Castor SKG. Two of them (RasTyc 0348+5202 
and RasTyc 2100+4530) have probability exceeding 80% which 
leaves little doubt as their membership to the Pleiades SKG. 
With probability higher than 50%, RasTyc 0307+3020 is an- 
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Fig. 10. The U-V (left panel) and U-W {right panel) kinematic diagrams for our moving group candidates (#i symbols ; see 
Table [6]for cross-match). The average velocity components (dots) of some young SKG and those of some late-type stars members 
of these young SKG are also plotted (square, triangle, circle, upside down triangle, and U symbols for the IC 2391 supercluster, 
Pleiades, Castor, UMa moving groups, and Hyades supercluster, respectively). The locus of the young-disc (YD) and old-disc (OD) 
populations are also marked. 



Table 6. Kinematics data and membership probability of moving group candidates. 



# Name 


RasTyc Name 


U 


v© 


w© 


IC 2391 


Pleiades 


Castor 


UMa 


Hyades 






(km s" 1 ) 


(km s- 1 ) 


(kms- 1 ) 












#1 BD+45 598 


RasTyc 0221+4600 


6.7 


-10.7 


-5.3 


N (0.0) 


Y(17.6) 


Y (65.4) 


N (0.0) 


N (0.0) 


#2 BD+29 525 


RasTyc 0307+3020 


12.3 


-19.6 


-17.0 


N (0.9) 


Y(54.1) 


N (0.9) 


N (0.0) 


N (0.0) 


#3 HD 275148 


RasTyc 0319+4212 


0.2 


-32.9 


-8.0 


N (0.0) 


N (9.8) 


N (0.0) 


N (0.0) 


N (0.0) 


#4 HD 22179 


RasTyc 0335+3113 


16.4 


-18.4 


-6.6 


Y (28.5) 


Y (48.5) 


N (6.5) 


N (0.0) 


N(1.2) 


#5 HD 23524 


RasTyc 0348+5202 


10.5 


-18.7 


-5.3 


N (0.6) 


Y(82.1) 


N (6.6) 


N (0.0) 


N (0.0) 


#6 HD 170527 


RasTyc 1825+6450 


40.4 


-58.3 


-5.9 


N (0.0) 


N (0.0) 


N (0.0) 


N (0.0) 


N (0.0) 


#7 BD+44 3670 


RasTyc 2100+4530 


11.0 


-22.8 


-8.2 


N (0.8) 


Y (99.3) 


N (0.6) 


N (0.0) 


N (0.0) 



other highly probable new Pleiades SKG member. The results 
are more ambiguous for RasTyc 0335+3113, which is either 
a new IC 2391 or Pleiades SKG member with probability 
close to 30 and 50 %, respectively. RasTyc 0221+4600 is our 
only candidate probably associated with the Castor SKG, al- 
though again we cannot rule out its connection to the Pleiades 
SKG. Unsurprisingly, the quantitative analysis confirmed that 
none of the stars classified as giants (RasTyc 0319+4212 and 
RasTyc 1825+6450) are associated with any of the SKG. 

The very good Hipparcos parallaxes (o- K /n < 0.1) of 
RasTyc 0348+5202 (#5) and RasTyc 1825+6450 (#6) allowed 
their positions to be plotted in the HR diagram, as illustrated in 
the right panel of Fig. [7] As expected, RasTyc 1825+6450 lies in 
the giant clump and is compatible with a 2.5 solar mass progen- 
itor ~ 500 Myr old according to post evolutionary tracks. On the 
other hand, RasTyc 0348+5202 is compatible with a ~ 30 Myr 
old solar mass PTTS evolving towards the ZAMS. These two 
particular cases for which we can access their evolutionary sta- 
tus, based on independent methods, give us a high confidence on 
our analysis and results. 



6. Conclusions 

In this paper we have presented the first results for a high- 
resolution spectroscopic survey of optical counterparts of 
X-ray sources in the northern hemisphere aiming at a thorough 
characterisation of the presumably youngest field stars in the 
solar neighbourhood. In particular, we focused on the >- 400 
optically brightest candidates for which we acquired about 800 
spectra in the Ha and/or lithium spectral regions. Using properly 
developed methods, we determined accurate radial and rota- 
tional velocities, effective temperature, and gravity. Moreover, 
an evaluation of the age (based on the Li 1/16707.8 line depth) 
and of the chromospheric activity level were performed. We 
have also identified several new binaries and multiple systems. 
We discussed the fraction of these systems in stellar soft X-ray 
surveys and showed that the contamination by lithium-rich 
evolved stars is not negligible. This sample of optically-bright 
and X-ray flux-limited sources turns out to mainly be composed 
of young stars (age < 1 Gyr) with a smaller contribution from 
an older population that was probably born within the last 
giga-years. Outstanding results are the discovery of about 50 
new lithium-rich giants and new members of young moving 
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groups, which are also PTTS candidates. 

This sample promises to give new insight into various fields 
of research, such as: i) the galactic disc heating and the recent 
local star formation rate, ii) the discovery of isolated very young 
stars that are promising targets for planet or debris-disc searches 
with ground-based or space telescopes, Hi) the evolution of 
binaries and multiple systems. 
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Appendix A: The data 

The data are published in electronic form only. The tables, avail- 
able at the CDS, contains all data for the 426 stars we observed 
at high resolution and discussed in this paper. The master tables 
(one for S, and SB1; the second one for SB2, and SB3) list all 
the time-independent data. In particular we have tabulated: the 
name (RasTyc and SIMBAD; columns #1 and #2), Tycho right 
ascension, and declination (columns #3 and #4), B - V colour 
index (column #5), ROSAT PSPC count rate, and associated er- 
ror (columns #6 and #7). In addition to these basic data, we also 
list the parameters derived in this study (but the radial veloc- 
ity): multiplicity status (column #8 : S for single or SB1, SB2, 
and SB3 for single-lined, double-lined, and triple systems, re- 
spectively), rotational velocity (column #9), effective tempera- 
ture (column #11), gravity (column #13), and metallicity (col- 
umn #15) with their associated errors (column #10, #12, #14, 
and #16, respectively). Finally we also list the measured net Ha 
(column #19 and #20) and lithium (column #21 and #22) equiv- 
alent widths and errors for single and SB1 stars only. For the 
case where the CCF had many peaks, we list the individual mea- 
surements in separate rows and refer to the primary (suffix a) as 
the component with the stronger absorption lines (i.e. the larger 
area of the CCF dip). The radial velocity measurements and er- 
rors are listed in dedicated tables together with the julian date of 
observation. 

Appendix B: Notes on PMS-like candidates 

We briefly summarize the known characteristics of these stars 
as recovered from the SIMBAD database and references therein. 
All the stars were observed by 2MASS. 

- RasTyc 0221+4600 : classified as a K0 star. 

- RasTyc 030 7+3020 : consid ered as a T Tauri candidate 
(WTTS) by [Li & Hul ([19981) o n the basis of the ROSAT 
RASS-BSC dVoges etal.lll996l) data selected by high hard- 
ness ratios. They quoted a spectral type G5IV. The star was 
also classified as PMS by iDucourant et al.l d2005l) in their 
PMS stars Proper Motion Catalogue. 

- RasTyc 0319+4212 : classified as a G5 T T auri star. 

- RasTyc 0335+3113 : classified as a GO star. iLi & Hul (fl998l) 
considered this star as WTTS candida te with spectral type 
G5IV. Also considered as PMS in the lDucourant et al.l d2005) 
catalogue. 

- RasTyc 0348+5202 : quoted as a double or multiple G6IV 
st ar dB artkevicius & Gudas] l2002f) while classified as KOV 
in lJacobv et al I dl984l) spectral atlas. ISilva & CorneiH d 19921) 
had taken a low resolution spectrum in the visible and near 
IR included in the New library of stellar op tical spectra. 

- RasTy c 1825+6450 : classified as G5IV bv lGunn & Strvker! 
d 1983b and identified as a ne w periodic variable from the 
Hipparcos epoch photometry dKoen & Everll2002l). I n their 
kinematic and dynamical studies. lFamaev et a l. (2005) found 
rich small-scale structure, with several clumps correspond- 
ing to various streams and superclusters. This star is part of 
their sample but it is not clear to what group they decided to 
attach it. 

- RasTyc 2100+4 530 : quoted as an F8 star. Observed at low 
resolution by Mot ch etall d 19971) in the framework of an 
identification programme of RASS galactic plane sources. 
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Fig. 11. Ha (Ze/f panel) and lithium (n'g/if panel) spectra (black solid line) of the new moving-groups PTTS candidates 
RasTyc 0221+4600, RasTyc 0307+3020, RasTyc 0335+3113, RasTyc 0348+5202, and RasTyc 2100+4530 (from top to bottom), 
together with the synthetic spectrum built up with the reference spectra broadened at the v sin i of the star's components and Doppler- 
shifted according to its RV (thin red line). 
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Fig. 12. Ha (left panel) and lithium (right panel) spectra (black solid line) of lithium-rich giants RasTyc 0319+4212 and 
RasTyc 1825+6450 (from top to bottom), together with the synthetic spectrum built up with the reference spectra broadened at 
the v sin i of the star's components and Doppler-shifted according to its RV (thin red line). 



